Griffiths-McCoy singularities in random quantum spin chains: 
Exact results through renormalization 
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The Ma-Dasgupta-Hu renormalization group (RG) scheme is used to study singular quantities in 
the Griffiths phase of random quantum spin chains. For the random transverse-field Ising spin chain 
we have extended Fisher's analytical solution to the off-critical region and calculated the dynamical 
exponent exactly. Concerning other random chains we argue by scaling considerations that the 
RG method generally becomes asymptotically exact for large times, both at the critical point and 
in the whole Griffiths phase. This statement is checked via numerical calculations on the random 
Heisenberg and quantum Potts models by the density matrix renormalization group method. 
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In a random quantum system at zero temperature sev- 
eral physical quantities are singular not only at the crit- 
ical point, but in a whole region, as well, which extends 
on both sides of the transition point [jj] . In this Griffiths 
phase |Q] a random quantum system is non-critical in the 
space direction (spatial correlations decay exponentially) , 
whereas it is critical in the time direction and the corre- 
sponding behavior due to Griffiths-McCoy singularities 
P,pJ is controlled by a line of semi-critical fixed points 
characterized by the dynamical exponent z(S), which de- 
pends on the value of the quantum control-parameter, 
5. For example average autocorrelations decay in (imag- 
inary) time, r, as ~ t~ z ^; for a small magnetic field, 
H — > 0, the magnetization behaves as ~ H 1 / 2 ^; the 
low temperature susceptibility and specific heat are sin- 
gular as ~ T 1 /^)" 1 and ~ T 1 / 2 ^, respectively. Let us 
mention that Griffiths-McCoy singularities are relevant 
in experiments on quantum spin glasses PJ and provides 
a theoretical explanation about the non-Fermi liquid be- 
havior in U and Ce intermetallics Q . 

Among the theoretical methods developed to study 
random quantum systems the renormalization group 
(RG) scheme introduced by Ma, Dasgupta and Hu 
Pj plays a special role. For a class of systems, the 
critical behavior of those is controlled by an infinite- 
randomness fixed point (IRFP), the RG method be- 
comes asymptotically exact during iteration. For some 
one-dimensional models, random transverse-field Ising 
model [| (RTIM) and the random XXZ-model @, Fisher 
has obtained analytical solution of the RG equations 
and in this way many new exact results and new phys- 
ical insight about the critical behavior of these models 
have been gained. Subsequent analytical Jl0[ and nu- 
merical jlf],[L0| investigations of the models are in agree- 
ment with Fisher's results. The RG-scheme has been 
numerically implemented in higher dimensions |^JT^|, as 
well, to study the critical behavior of the RTIM and rea- 



sonable agreement with the results of quantum Monte- 
Carlo simulations |f3) has been found. Considering the 
Griffiths-phase of random quantum spin chains here the 
RG-scheme has been rarely used p^ |, mainly due to the 
general belief that the method looses its asymptotically 
exact properties by leaving the vicinity of the scale in- 
variant critical point. 

Our aim in the present Letter is to clarify the appli- 
cability of the Ma-Dasgupta-Hu RG-method in the Grif- 
fiths phase of random quantum spin chains. We start 
with the RTIM, extend Fisher's calculation || to the 
Griffiths phase of the model and present the analytical 
solution of the RG-equations. Then, for general mod- 
els, we analyze by scaling considerations the structure of 
the RG equations around the line of semi-critical fixed 
points and arrive to the conclusion that the RG method 
becomes asymptotically exact in the whole Griffiths re- 
gion. This statement is then checked numerically on the 
random Heisenberg chain and the random quantum Potts 
model (RQPM) using the density matrix renormalization 
group (DMRG) method. 

We start with the Id RTIM which is defined by the 
Hamiltonian: 



Hi = - -^X- i - 53 hi<T i 



(1) 



where the af' z are Pauli matrices at site i and Ji and hi 
are the couplings and the transverse fields, respectively, 
which are independent random variables. The quantum 
control-parameter of the model is defined as 



[ln/i] av -[lnJ] 

av 

var[ln h] + var[ln J] 



(2) 



where var[a;] is the variance of x and we use [. . .] av to de- 
note averaging over quenched disorder. For S > (S < 0) 
the system is in the paramagnetic (ferromagnetic) phase, 
so that the random quantum critical point is at 5 = 0. 
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In the Ma-Dasgupta-Hu RG-method the strongest 
term in the Hamiltonian, coupling or field, of strength 
Cl are successively decimated out and the neighboring 
fields or couplings are replaced by weaker ones, which 
are generated by a perturbation calculation. The ba- 
sic RG-equations for coupling and field decimations are 
given by: 



hih 



i'H+l 



J 



Jj-\Jj 

hjK 



(3) 



respectively, which are related through duality. Here, for 
the RTIM k = 1. Under renormalization we follow the 
probability distributions of the couplings, R(J, Cl), and 
that of the fields, P(h, Cl). When the energy-scale is low- 
ered as Cl — * Cl — dCl the distribution of the couplings is 
changed as: 



dR(J, H) 

dn 



= r(j, n)[p(n, n) - r{ci, ci)} 



. JCln „, CIk 



-p(ci,ci) dfR(j',n)R(^,ny-^ , (4) 

where the first term in the r.h.s. is due to a balance be- 
tween decimated couplings and normalization, whereas 
the second term accounts for the generated new cou- 
plings. A similar equation for the field distribution fol- 
lows from Eq.(^|) through duality 



dP{h,Cl) 



an 



= p(h, ci)[r(ci, ci) - p(ci, ci)} 



hCln „, CIk 



r(ci,ci) dh'p(h',n)P(^,n)-£ 



(5) 



which amounts to interchange J <-> ft, and i? <-> P. 

For the RTIM, i.e. for /t = 1, we found one class of 
solutions of the RG equations in the form: 



R(J, Cl) = R(Cl, CI) (Cl/J) 1 - 11 ^'^ 

p{h, ci) = p(ci, ci) (n/h) 1 ~ p(p,a)a 



(6) 



where the distributions involve the parameters, R(Cl) = 
R(C1,CI) and P(Ct) = P(Cl,CY), which satisfy the relation 
(P — R)C1 = l/z = const. Thus the solution is charac- 
terized by one parameter, z — z{5), which depends on 
the quantum control-parameter, 5: at the critical point, 
1/2(0) = 0, whereas in the paramagnetic phase, S > 0, 
l/z(S) > and monotonically increases with S. 

In terms of the variables, y = RCt + 1 /2z = PCI — 1 jlz 
and x = — In CI we obtain the differential equation: 



dy 
dx 



y 



1 

4^2 



(7) 



which has the solution, y — \/(x — Xq), xq = const, at the 
critical point with l/z = 0. The distribution of p = RCl 
in terms of the variable 77 = — (In CI — In J)/ In CI is given 



by p(?7)d?7 = exp(— n)dn, which corresponds to the fixed- 
point solution by Fisher j8j. At this point we refer to 
Fisher's analysis || and conclude that the functions in 
Eqs(^) indeed represent the fixed-point distribution for 
all non-singular initial distributions. 

In the Griffiths phase, S > 0, the solution of Eq.(^J) in 
terms of the original energy-scale variable, Cl, is given by 



y 



y /2z + 1 /4z 2 tanh [\n{Cl /Cl) /2z] 
1 /2z + y tanh [ln(Cl /Cl) ~/2z) 



(8) 



where y = yo at a reference point, Cl — CIq- Approaching 
the line of semi-critical fixed points, i.e. for CI/CIq — > 0, 
we have to leading order: 



RCl = PCI- l/z = R(Cl )/[P(Cl )z] (O/Oq 



■V* 



(9) 



thus P and R have different low energy asymptotics. 

The physical relevance of l/z can be obtained by study- 
ing the change of number of spins, n^ — > nji — dn^ 
connected with a change in the energy scale as Cl — > 
Cl — dCl. This leads to the differential equation dnn/dCl = 
Tin [P(Cl, Cl) + R(Cl, Cl)), the solution of which is given by: 

nn = {cosh[ln(fi /fi)/2z] + 2z y sinh [ln(Cl /Cl)/2z]y 2 , 

(10) 

which along the line of semi-critical fixed points has 
the asymptotic behavior nn — const Cl 1 /*, Cl — > 0. 
Since the typical distance between remaining spins is 
Ln ~ 1/fiQ ~ Cl~ x / Z , we can identify z as the dynam- 
ical exponent, which governs the relation between time- 
and length-scales as t ~ L z . 

Next we show that z is invariant along the RG trajec- 
tory and can be deduced from the original distributions. 
For this we consider the averages, [J M ] av and [h~^] &v , and 
using Eqs.(0) and (0) we calculate the derivative: 



k-wl = a - imru 

+ R(Cl,Cl)Cl»[h-»] 



dCl 

x (P(Cl,Cl)Cl-v{J. 



(11) 



which is vanishing for /1 — jx, if [(J//i)' i ] av = 1. Conse- 
quently fl stays invariant along the RG trajectory until 
the fixed point, where using the distribution in Eqs.(|^) 
we obtain fl — l/z. Thus the dynamical exponent for the 
RTIM is given by the solution of the equation: 



(J/h) 



= 1 



(12) 



which is then exact, since the RG-transformation be- 
comes asymptotically exact as Cl — > 0. This latter state- 
ment follows from the fact that the ratio of decimated 
bonds, Anj, and decimated fields, An/ l; goes to zero as 
Anj/An h = R{Cl,Cl)/ P(CI,C1) - Cl 1 / 2 . Then the proba- 
bility, Pr(a), that the value of a coupling, J, being neigh- 
bor to a decimated field is Cl < J < aCl with < a < 1 
is given by Pr(a) — 1 — a , which goes to zero for any 
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non-zero a, since Pf2 ~ &}/ z at the fixed point. Con- 
sequently the decimations in Eq.(||) and the related RG 
equation in Eqs.(0) and (@) are indeed exact. 

The relation in Eq. (|l^) can be recovered starting from 
the exact expression for the surface magnetization, m s , 
on a chain with a fixed spin at site L+l JlJ|,[l(| : 



En 

1=1 i=l 



(13) 



This type of expression is the so-called Kesten-variable 
in the mathematical literature |H| and the correspond- 
ing probability distribution is singular at m s = in the 
thermodynamic limit for [In J] av > [ln/i] av [^5|,^G|. Then 
for the distribution of In m s we have the following singu- 
larity: 



P(lnm s )~my z , m s — > 



(14) 



where z is the solution of Eg. (|l2j) in terms of the dual 
variables. The physical origine of the small m s tail of the 
distribution is due to such samples which have a weakly 
coupled domain (WCD), which effectively cuts the sys- 
tem inty two very weakly interacting parts and thus re- 
duces the surface order enormously. In the dual system 
in the paramagnetic phase the dual object to a WCD is 
a strongly coupled domain (SCD) which results in a very 
small energy gap, e. Thus in the tails of the distributions, 
m s and e are dual quantities. This remark indicates that 
the dynamical exponent found by the RG calculation in 
Eq.Q is indeed exact JL7|. 

Next, we consider general random quantum spin chains 
with a critical IRFP and analyze the structure of the RG 
equations close to the line of semi-critical fixed points, 
thus as il — > 0. As for the RTIM, the decimation for fields 
and couplings is asymmetric and for O — > exclusively 
fields are decimated out, which are typically infinitly 
stronger, than the neighboring couplings. Therefore the 
RG decimation equations in Eq.([|) are asymptotically 
exact. The second point is to show that the dynamical 
exponent stays invariant along the RG trajectory, even 
though in the starting phase the RG equations are ap- 
proximative. For this we consider the low energy tail of 
the distribution function of the first gap, Pi(lne), which 
involves in analogy to Eq. ( 14 ) the exponent z, and use the 
scaling result of Rcf |l^] . This states that the probability 
distribution of the second, third, etc. gaps are related 
to Pi (Lie) as P 2 (lne) - Pi(hie), P 3 (hie) ~ if (me), 
etc, due to the fact that for a small second, third gap 
one needs two, three independent SCD-s and the corre- 
sponding probabilities are multiplied. In the RG deci- 
mation the SCD-s are only eliminated through coupling 
decimation, since their couplings are stronger than the 
average fields. If at some time a SCD with a small 
gap, e, is eliminated then in the probability distribu- 
tion, Pi(lne), one should consider the former second 



gap and use the corresponding conditional probability, 
Pi (Lie) -> P 2 (lne)/Pi(lne) - Pi (Lie). Thus the small 
energy tail of the gap-distribution and consequently the 
dynamical exponent remains invariant under the renor- 
malization procedure. The previously obtained exact re- 
sults for the RTIM give strong support for the validity of 
these phenomenological considerations. 

For a numerical demonstration of the validity of the 
above statement we considered two random quantum 
spin chains, the dimerized Heisenberg (XXX) chain and 
the q-state RQPM, both having a set of RG equations 
very similar to that of the RTIM in Eq.(|). For the 
dimerized XXX chain J and h in Eq.(|^) are replaced by 
the Heisenberg couplings at odd and even positions, J Q 
and J e , respectively, and the parameter takes the value 
re = 2 HJ5). The distance from the critical point is mea- 
sured similarly to Eq.(0). For the q-state RQPM fields 
and couplings play analogous role as for the RTIM, the 
quantum control-parameter is given in Eq.(^J), whereas 
re takes the value k = q/2 |l9| . We note that the RG 
equations for the XXX-chain and the q = 4 state RQPM 
are identical. 
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FIG. 1. Dynamical exponents from numerical iteration of 
the RG-equations in Eq.(H) for different values of the param- 
eter re. For k = 1, which corresponds to the RTIM and the 
random dimerized XX-chain the exact result is given by the 
full line, for re — 2 and 4 the broken lines are guide to the eye. 

At the critical point the RG-equations for 1 < K < oo 
has been solved by Senthil and Majumdar |l^| with the 
result that re is an irrelevant variable and the IRFP is the 
same as for the RTIM. In the Griffiths-phase we could not 
find a complete solution of the RG equations, in spite of 
the close similarity to that of the RTIM. We could, how- 
ever, show that up to an accuracy of 0{VL x / z ) the solution 
is of the form of Eqs.(Q) and thus there is infinite ran- 
domness along the line of fixed points. The z exponent, 
however, does depend on the parameter re, since the va- 
lidity of the condition in Eq. (O) is limited to re = 1 , thus 
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in general z = z K [S). 

We have calculated the dynamical exponent by a nu- 
merical implementation of the RG scheme over 50000 
samples of length L < 2 14 . Starting with the uni- 
form probability distribution: Rq(J) = 0(1 — J)Q(J), 
Po(h) = Q{ho — h)Q(h)/h (and analogously for the 
XXX-chain) we got the estimates shown in Fig.l: l/z K is 
a monotonously decreasing function of k and eventuclly 
it goes to zero in the whole Griffiths phase in the limit 




FIG. 2. Probability distribution of the logarithm of the 
surface magnetization of the dimerized random XXX-chain 
with S = — (ln2)/2 for different finite systems calculated 
by the DMRG method. The full straigth line with a slope 
l/z(5) = .47 represents the asymptotic behavior according to 
the RG method, as given in Fig. 1. In the inset the same 
quantity is shown for the q = 4 state RQPM at the same 
distance from the critical point. The asymptotic slope of this 
distribution can be well fitted by the same exponent as for 
the XXX-model. 

The dynamical exponents of the XXX-chain and the 
RQPM are also calculated directly from the asymptotic 
behavior of the distribution of the surface magnetization, 
as given in Eq.(p^). For the numerical calculations of the 
surface magnetizations we used the DMRG method for 
rather large finite chains with L < 64 and considered 
some 20000 samples. We found an overall agreement be- 
tween the dynamical exponents calculated by the two 
methods. As a demonstration we show in Fig. 2 the dis- 
tribution of m s for the XXX-chain, compared with that 
of the q = 4 state RQPM, where for both models we are 
at the same distance from the transition point. As seen 
in Fig. 2 the asymtotic behavior of the two distributions 
is identical, as expected on the RG basis, since k = 2 
for both models. Furthermore the dynamical exponents 
agree very well with those calculated by the RG method. 

To conclude, we have shown in this paper that for ran- 



dom quantum spin chains the RG method of Ma, Das- 
gupta and Hu is asymptotically exact for large times, 
i.e. along the line of semi-critical fixed points. Conse- 
quently the dynamical exponent in the Griffiths phase 
calculated by the RG-method is exact. This result will 
hopefully stimulate new efforts to solve the RG equations 
for different systems analytically in Id also outside the 
critical point, and in higher dimensions to perform nu- 
merical calculations and get more precise estimates for 
the Griffiths-McCoy singularities. 
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